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The antiferromagnetic Fulde-Ferrell-Larkin-Ovchinnikov (AFM-FFLO) state of coexisting 
d-wave FFLO superconductivity and incommensurate AFM order is studied on the basis of 
Bogoliubov-de Gennes (BdG) equations. We show that the incommensurate AFM order is 
stabilized in the FFLO state by the appearance of the Andreev bound state localized around 
the zeros of the FFLO order parameter. The AFM-FFLO state is further enhanced by the 
induced 7r-triplet superconductivity (pair density wave). The AFM order occurs in the FFLO 
state even when it is neither stable in the normal state nor in the BCS state. The order 
parameters of the AFM order, d-wave superconductivity, and 7r-triplet pairing are investigated 
by focusing on their spatial structures. Roles of the spin fluctuations beyond the BdG equations 
are discussed. Their relevance to the high-field superconducting phase of CeCoIn 5 is discussed. 

KEYWORDS: FFLO state, antiferromagnetism, 7r-triplet pairing, CeColn 5 



1. Introduction 

The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 
in superconductors was predicted in the 1960s by Fulde 
and Ferrell, 1 and Larkin and Ovchinnikov. 2 In contrast to 
the Bardeen-Cooper-Schrieffer (BCS) state, Cooper pairs 
have a finite total momentum in the FFLO state, which 
leads to the spontaneous breaking of the spatial sym- 
metry. Although this novel superconducting state with 
an exotic symmetry has been attracting much interest, 
the experimental search for this state had been fruit- 
less for nearly 40 years. Under these circumstances, the 
discovery of a new high-field superconducting (HFSC) 
phase in CeCoIns, 3,4 which is a likely candidate for the 
FFLO state, triggered many theoretical and experimen- 
tal studies. 5 This recent interest on the FFLO supercon- 
ductivity/superfluidity extends further in various related 
fields, such as organic superconductors, 6-10 cold atom 
gases, 11,12 astrophysics, and nuclear physics. 13 

The HFSC phase of CeCoIns has been interpreted 
widely within the concept of the FFLO state. 5 ' 14-24 How- 
ever, recent observations of the AFM order in the HFSC 
phase call for a reexamination of this conclusion. 25-27 
It is expected that this AFM order will be closely re- 
lated to the AFM quantum critical point observed in 
CeCoIns. 28,29 Therefore, the nuclear magnetic resonance 
(NMR) 25,27 and neutron scattering 26 measurements may 
have uncovered a novel superconducting state in this 
strongly correlated electron system. In particular, the 
neutron scattering measurement has explored the prop- 
erties of the AFM order and found that the wave vector 
of the AFM order is incommensurate q\c = Q+5\q with 
8\c along the [1,-1,0] direction for magnetic fields in the 
a6-plane of the tetragonal lattice. 26 The AFM staggered 
moment Maf is directed to the c-axis. 

Some scenarios have been proposed for the HFSC 
phase of CeCoIn 5 from the theoretical point of view. We 
have investigated the AFM order arising from the nodal 



quasiparticles in the inhomogeneous Larkin-Ovchinnikov 
state. 30 The SDW order triggered by the emergence 
of 7r-triplet pairing 31-34 has been investigated in the 
BCS state 31,32,34 and in the homogeneous Fulde-Ferrell 
state. 33 In this study, we examine the AFM-FFLO state, 
in which the AFM order appears in the inhomogeneous 
FFLO state, on the basis of BdG equations. The typical 
phase diagram in the H-T plane and the spatial structure 
of the AFM-FFLO state are investigated in detail. 

2. Formulation 

Our theoretical analysis is based on the microscopic 
model, 

H=-t V ct c , +*' V ct c , 

Z ✓ i ,<j 3 i* 7 Z ✓ i ,a 3 ■><* 

<i ,j >,cr «i ,j >>,cr 

i <i,j> 

+ J E S r Sj -<7b#£%, (1) 

<i,3> i 

where is the spin operator and is the number op- 
erator at site i. To describe the quasi- two-dimensional 
electronic structure of CeCoIn 5 , we assume a square lat- 
tice, in which the bracket < i,j > (<< i,j >>) 
denotes the summation over the nearest-neighbor sites 
(next- nearest-neighbor sites) . The on-site repulsive inter- 
action is given by U, and V and J stand for the attrac- 
tive and AFM exchange interactions, respectively, be- 
tween nearest-neighbor sites. Wc assume V to stabilize 
the c?-wave superconducting state within the mean field 
BdG equations and choose J > for the AFM correla- 
tion in CeCoIns. We assume U — 1.15, V = —0.4, and 
J = 0.53 throughout this paper. The (i-wave supercon- 
ductivity and the significant AFM spin fluctuation can 
be self-consistently described using many-body theories 
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such as the fluctuation exchange (FLEX) approximation 
for the Hubbard or periodical Anderson model. 35 ' 36 We 
here assume the interactions V and J to describe these 
features in the inhomogeneous LO phase on the basis 
of BdG equations. Although the BdG equations neglect 
the AFM spin fluctuation beyond the mean field approx- 
imation, they are suitable for studying the qualitative 
features of the inhomogeneous superconducting and/or 
magnetic state. We will discuss the roles of AFM spin 
fluctuation at the end of this paper. 

With the last term in cq. (1), we include the Zeeman 
coupling due to the applied magnetic field. The g- factor 
is assumed to be g& = 2. We assume the magnetic field 
to lie parallel to the [100]-axis, which we choose to be 
the quantized axis of spin, with H — Hz. In this paper, 
we identify the a-, b-, and c-axes of the tetragonal lattice 
with the z-, x-, and ?/-axes for the spin, respectively. We 
choose the unit of energy such that t = 1 and t' jt — 
0.25. The chemical potential enters as \i — [i + {\U + 
AV)uq, where n is the number density for U = V = 
J = H = 0. We choose ^ = —1-05 so as to reproduce 
the incommensurate AFM order observed in the neutron 
scattering measurement for CeCoIn5. 26 Then, we obtain 
the number density n ~ 0.77. 

The BdG equations are formulated in a standard man- 
ner. We take into account the Hartree term arising from 
U, V, and J. The BdG equations are self-consistently 
solved for the mean fields of the spin < >, charge 
< n-r >, and superconductivity A% a -r =< , >. 

Since the magnetic moment is nearly opposite between 
the nearest-neighbor sites, the order parameter of the 
AFM order is described by the staggered moment de- 
fined as 

M A¥ {i) = (-1)™+" < S t >, (2) 

with % = (m,n). Two components of the pairing field 
A?°j play dominant roles in the following results. The 
first component is the c?-wave spin singlet pairing whose 
order parameter is described as 

A d (^) = AlK „ + At^ _-At^ .-Al^ (3) 

l,l+a a i ,i +b i ,i —b 

with a and b being the unit vectors along the a- and 
&-axes, respectively. The second component is the equal 
spin 7r-triplet pairing whose order parameter is described 
by the generalized d- vector d a ,b(0 as 

Af- ±s = ±(-l)™+^(-<7<(?) (4) 

A?* ?±S = ±(-l) m + nl -(-adl(I) + «£(?)). (5) 

The (i-wave spin singlet pairing is finite in the entire su- 
perconducting region, while the 7r-triplet pairing appears 
in the state of coexisting AFM order and superconductiv- 
ity. 3 1-33 ' 37-40 Because of the linear coupling between the 
magnetic, spin singlet pairing, and 7r-triplet pairing or- 
der parameters, two finite order parameters among them 
induce the other order parameter. Our model eq. (1) in- 
cludes the interactions leading to the magnetic order and 
spin singlet d-wave superconductivity, while the pairing 
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interaction for the 7r-triplet pairing is negligible. There- 
fore, the 7r-triplet pairing does not belong to the domi- 
nant orders, but is induced as a secondary order param- 
eter by the AFM order. 

The mean field Hamiltonian is obtained as 

H = -t V ct c la +t' V ct c , 

<i ,j >,cr «i ,j >>,a 

i $ 

+ £(-2tf < s, > +JJ2 < S US > -9bH)S~ 

i 5 

+ - V \A%< -ct ct ? +C.C.1, (6) 

where the summation of 8 is taken over 8 = ±a, ±6. The 
pairing field is obtained as A^,. - = — (V — J/4) < 

C7 tC? _l x i > —J/2 < Cj- iCv ? > and AZ 17 ^ ? = 

i ,J i+o,l 1 1 ,1 2+o,T 1,1+8 

— {V + J/4) < a c^ + g a >. The thermodynamic aver- 
age <> is obtained on the basis of the mean field Hamil- 
tonian eq. (6). The self-consistent equations for < >, 
< >, and A-°- - for each i yield the BdG equations. 
The free energy is obtained as 

F = - \t Iog[l + exp(-£7 a /T)] 

a 

+ J2(\ U < n i > 2 +l v H < n i >< n ns >) 

i 6 
+ J2(-U < S;> 2 + 1 -jJ2<S;>-< § US >) 

i s 

+1 y a%k <d - ci >, (7) 

2 i,t+8 i +6, a' i, a v ' 

where E a is the energy of Bogoliubov quasiparticles ob- 
tained by the mean field Hamiltonian eq. (6). The stable 
phase is determined by minimizing the free energy for 
the self-consistent solutions of the BdG equations. 

3. Results 

3.1 Phase diagram 

We first show the phase diagram against the temper- 
ature and magnetic field. The normal, uniform BCS, 
purely FFLO, and AFM-FFLO states are shown in 
Fig. 1. At high magnetic fields, the phase transition from 
the normal state to a superconducting state is of the 
first order, consistent with the experimental results for 
CeCoIns. 5,41,42 We have shown that both the on-site re- 
pulsion U and AFM interaction J are required to re- 
produce the first order phase transition to the FFLO 
state. 43 This indicates that both the AFM spin fluctua- 
tion and the local electron correlation play an essential 
role in the phase diagram of CeCoIn 5 . These features are 
qualitatively understood on the basis of the Fermi liquid 
theory. 44 Although the size of the numerical calculation 
(40 x 40 lattices) is not sufficiently large for distinguish- 
ing the first and second order phase transitions from the 
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BCS state to the FFLO state (BCS-FFLO transition), 
the BCS-FFLO transition is expected to be of the second 
order, as shown for a similar model. 30 The second order 
BCS-FFLO transition is described by the nucleation of 
the FFLO nodal plane. 30 ' 44 ' 45 Because of our tractable 
system size limitation, maximally two nodal planes fit 
into our calculation, as shown in Fig. 2. 

An important finding obtained from Fig. 1 is the ap- 
pearance of the AFM-FFLO state. The Neel temperature 
Tn in the BCS and normal states is less than 0.01, which 
is the lower limit of our calculation, while we obtain a 
much higher Neel temperature Tn of ~ 0.02 in the FFLO 
state. Thus, the AFM order is favored in the FFLO state 
rather than in the BCS and normal states. 
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Fig. 1. (Color online) Phase diagram for the magnetic field and 
temperature. The normal, uniform BCS, FFLO, and AFM- 
FFLO states are shown in the figure. Solid lines show the phase 
boundary between these states. The thin dashed line shows the 
fictitious phase transition line from the FFLO state to the AFM- 
FFLO state where the order parameter of 7r-triplet pairing is 
neglected. Note that the T c of superconductivity at H = is 
T c = 0.096. 



Two mechanisms stabilize the AFM order in the FFLO 
state. The first mechanism is the appearance of Andreev 
bound states around the spatial nodes of the modulated 
superconducting order parameter in the FFLO state. 
The 7r-phase shift of the order parameter introduced in 
the FFLO state produces a large local density of states 
(DOS) and triggers the AFM order. 30 This is the reason 
why the AFM order is favored in the FFLO state rather 
than in the normal and BCS states. The appearance of 
the Andreev bound states is a characteristic feature of 
the inhomogeneous Larkin-Ovchinnikov state and is ab- 
sent in the homogeneous Fulde-Ferrell state. The role of 
the Andreev bound states in the FFLO nodal structure 
is much more pronounced than that of the single vortex 
that has been investigated by Ogata. 46 This is because 
the FFLO nodal plane is a two-dimensional object, while 
the vortex is an one-dimensional object. 

The second mechanism is the linear coupling between 
the AFM order and the 7r-triplet pairing that favors the 
AFM order in the spin singlet superconducting state. It 
has been shown that this coupling stabilizes the AFM 
state in the BCS state 31 ' 32 and that in the Fulde-Ferrell 
state. 33 In contrast to the previous studies, 31 ' 32 the ir- 



triplet pairing state without the AFM order is hardly 
stabilized in our model. However, the coupling to the 
7r-triplet pairing significantly stabilizes the AFM order 
in the Larkin-Ovchinnikov state. The thin dashed line in 
Fig. 1 shows the fictitious Neel temperature in the FFLO 
state, if the 7r-triplet pairing is neglected. We observe the 
substantial increase in Tn due to the admixed 7r-triplct 
pairing. Thus, the 7r-triplet pairing plays a quantitatively 
important role in the stability of the AFM-FFLO state 
even when the amplitude of 7r-triplet pairing is small. 

We now turn to the question why the AFM order ap- 
pears discontinuously at the normal-FFLO transition as 
well as at the BCS-FFLO transition. The former discon- 
tinuity is simply caused by the first order phase transi- 
tion to the FFLO state. The order parameter of d-wave 
superconductivity appears discontinuous at the normal- 
FFLO transition and the AFM state is directly coupled 
to this order. The discontinuity at the BCS-FFLO tran- 
sition looks more surprising because the phase transition 
is of the second order. The key lies in the appearance 
of the Andreev bound states that we have mentioned 
above. The second order phase transition from the BCS 
state to the FFLO state is associated with the nucleation 
of domain walls in the superconducting order parame- 
ter. 30 ' 44,45 The Andreev bound states localized around 
such domain walls are the source of the AFM instabil- 
ity. Since the Andreev bound states in the isolated do- 
main wall arc nearly independent of the adjacent domain 
walls, within the mean field approach, the Neel temper- 
ature is nearly independent of the density of domain 
walls. In other words, the AFM order occurs immedi- 
ately when domain walls are nucleated at the BCS-FFLO 
transition. 30 Note that the spatially averaged magnetic 
moment is continuous at the second order BCS-FFLO 
transition, although the Neel temperature is discontin- 
uous there. The phase boundary near the BCS-FFLO 
transition obtained here by the BdG equations would be 
altered by the spin fluctuation, as discussed in §4. How- 
ever, the spin fluctuations do not alter the result that the 
AFM order is confined in the FFLO state. This feature 
of the phase diagram is consistent with the experiment of 
CeCoIn 5 that revealed the magnetic order in the HFSC 
phase, but neither in the normal state nor in the low field 
superconducting phase. 26 

3.2 Spatial structure of AFM-FFLO state 

We here investigate the spatial structures of the AFM- 
FFLO state shown in Fig. 1. Figure 2(a) depicts the order 
parameter of the d-wave spin singlet superconductivity 
A d (i ), while Fig. 2(b) shows the AFM staggered moment 
M^ F (i ) for i = (m, n). In Fig. 2(a), we assume that the 
modulation vector of FFLO superconductivity is paral- 
lel to the magnetic field q fflo II H, in accordance with 
ref. 22. We find that the AFM staggered moment is per- 
pendicular to the applied magnetic field, consistent with 
the neutron scattering measurement. 26 Since we neglect 
the spin-orbit coupling, the AFM state with Maf(j ) || x 
is degenerate with that with Maf(0 II V- We choose the 
former in the following results. 

We observe that the AFM staggered moment is local- 
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Fig. 2. (Color online) Spatial dependence of (a) the spin singlet 
pairing field with d-wave symmetry A d (i) and (b) the AFM 
staggered moment M^ F (i ) for i = (m,n). We assume T = 0.01 
and H = 0.112. The direction of the magnetic field is shown by 
an arrow. 



ized around the spatial nodes of the modulated order 
parameter of d-wave superconductivity. This is because 
the Andreev bound states mainly induce the AFM or- 
der, as previously mentioned. Another intriguing find- 
ing is the direction of the incommensurate wave vector, 
which is perpendicular to the FFLO modulation vector 
as 5\c -L q fflo II H. This structure is consistent with 
the experimental result of CeCoIns. 26 The amplitude of 
the incommensurate wave vector | <5 io I is independent of 
the density of FFLO nodal planes, as we have shown in 
ref. 30. Therefore, our results are compatible with the 
experimental observation in which q iq = Q + 8 ic is in- 
dependent of magnetic field. 26 Thus, our results on the 
direction, spatial structure, and magnetic field depen- 
dences of the AFM staggered moment are consistent with 
the neutron scattering measurement. 26 

Figure 3 shows the order parameter of the 7r-triplet 
pairing. We obtain the generalized d-vector as d a (i) ~ 
— G?b(0 oc x + iay with < a < 1. This structure 
is the same as the d-vector proposed for the high-field 
superconducting phase in Sr2Ru04. 47 Since d y ah (i) is 
pure imaginary for any i , we plot the real quantity 
— d*(i) + id^(i) in Fig. 3. The complex spatial struc- 
ture of the 7r-triplet pairing arises from the spatial mod- 
ulation in the AFM staggered moment M£ F (z) along 
the 6-direction and that in the spin singlet pairing field 
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Fig. 3. (Color online) Spatial dependence of the order param- 
eter for the 7r-triplet pairing. We plot — dj(i) + ida(i) = 
- At T , _). The directions of the FFLO mod- 
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ulation q fflo an d incommensurability S ic are shown by ar- 
rows. The parameters are the same as those in Fig. 2. 



A d (i ) along the a-direction. Since the 7r-triplet pairing is 
induced by the combination of the AFM order and spin 
singlet pairing, the generalized d- vector for the 7r-triplet 
pairing d a ,b(0 changes the sign at the zeros of M^ F (i) 
and those of A d (z ). The modulation wave vectors of the 
7r-triplet pairing are mainly q fflo along the a-axis and 
Sic along the &-axis. 
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Fig. 4. (Color online) Spatial dependence of the magnetization 
M z (i) =< S% > parallel to the magnetic field. The parameters 
are the same as those in Fig. 2. 



Finally, we show the spatially inhomogeneous magne- 
tization along the magnetic field M z (i) (Fig. 4). The 
magnetization mainly arises from the polarization of the 
Andreev bound states. 48 Therefore, we clearly observe 
the magnetization localized around the spatial nodes of 
the spin singlet pairing field. In addition, a weak mod- 
ulation of magnetization appears along the nodal lines 
owing to the incommensurate AFM order. It is shown 
that the magnetization is enhanced at the intersection 
points between the nodal planes of the incommensurate 
AFM order and the spin singlet pairing. 
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3.3 Neutron scattering 

We here propose an experiment that can unambigu- 
ously identify the AFM-FFLO state. Figure 5 shows the 
Fourier transformation of the magnetic moment perpen- 
dicular to the magnetic field, 

S x (q) = ^Y / <St>eMW-*)- (8) 

i 

We observe two sharp peaks in S x (q) at the incommen- 
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q x /2n 

Fig. 5. (Color online) S x (<?) in the AFM-FFLO state for q 
(q x ,q y ). The parameters are the same as in Fig. 2. 



surate wave vectors qic = Q ± 5\c where 8 ic /2n = 
(0,0.075) in our calculation. In addition to these main 
peaks, the satellite peaks appear at q = qic + <f s t with 
q st || H _L Sic in the AFM-FFLO state. These satellite 
peaks are the characteristic feature of the AFM-FFLO 
state and unambiguously show the spontaneous transla- 
tion symmetry breaking along the magnetic field. There- 
fore, if the spin structure in Fig. 5 is observed by a neu- 
tron scattering measurement, that would give unambigu- 
ous evidence for the AFM order in the FFLO state. The 
neutron scattering measurements have revealed a sharp 
peak at q = qic, 26 however, the presence of satellite 
peaks shown in Fig. 5 has not yet been explored experi- 
mentally 

4. Roles of Spin Fluctuation 
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Fig. 6. (Color online) (a) Schematic view of the H-T-pha.se di- 
agram obtained by the mean field BdG equation in this paper, 
(b) Phase diagram in which the spin fluctuation is phcnomeno- 
logically taken into account. The solid line in the superconduct- 
ing state shows the BCS-FFLO transition, while the dashed line 
shows the Neel temperature. 



We turn to the roles of the spin fluctuations ne- 
glected in the BdG equations. First, the curvature of 
the BCS-FFLO transition line is changed by spin fluc- 
tuations near the AFM quantum critical point. We have 
investigated the FFLO superconductivity in the two- 
dimensional Hubbard model on the basis of the FLEX 
approximation. The resulting BCS-FFLO transition line 
shows a concave curvature, 36 as observed for the HFSC 
phase of CeCoIns. 19 This is mainly due to the quasipar- 
ticle lifetime renormalization by the spin fluctuation that 
suppresses the FFLO state in the high-temperature re- 
gion. 36 This is in sharp contrast to the BdG equation 
that shows the convex BCS-FFLO transition line as in 
Fig. 1. 

Second, we discuss the mechanism that stabilizes the 
incommensurate AFM order in the FFLO state. Spin sus- 
ceptibility in the spatially uniform state is enhanced at 
q — Q + 5 ic nearly independent of the direction of 5 ic 
in the ah- plane on the tetragonal lattice. Note that the 
directional fluctuations of 5 ic suppress the weakly in- 
commensurate AFM order in the uniform phases, such 
as the normal and BCS states. On the other hand, the 
appearance of the significant in-plane anisotropy due to 
the modulated FFLO order parameter suppresses the di- 
rectional fluctuations and favors the AFM order. This 
mechanism may play a role in stabilizing the AFM order 
in the HFSC phase of CeCoIns . 

Finally, the spin fluctuations play another role just 
above the BCS-FFLO transition line. The continuous 
phase transition from the BCS state to the FFLO state 
is described by the nucleation of domain walls. The den- 
sity of domain walls approaches zero near the BCS- 
FFLO transition line. Then, the spatial dimensionality 
of the AFM order is reduced from three to quasi-two- 
dimensions. The reduced dimension enhances the fluctu- 
ations and suppresses the long-range order at finite tem- 
peratures when we neglect the broken SU(2)symmctry 
due to the spin-orbit coupling. The AFM order is sup- 
pressed by this effect just above the BCS-FFLO transi- 
tion line. The density of domain walls rapidly increases 
with growing magnetic field from the BCS-FFLO transi- 
tion line, and then, this effect of spin fluctuation is sup- 
pressed. 

Taking into account the roles of spin fluctuations, we 
obtain the schematic phase diagram in Fig. 6(b). The 
AFM order is confined in the FFLO state, but it is sup- 
pressed around the BCS-FFLO transition line. The pure 
FFLO state is stabilized just above the BCS-FFLO tran- 
sition line. The two phase transition lines, namely, the 
BCS-FFLO transition and the AFM order, are close to 
each other. Therefore, it may be difficult to distinguish 
these transition lines in an experiment. The CeCoIn 5 at 
ambient pressure seems to be the case. 

5. Summary and Discussion 

We investigated the incommensurate AFM order in the 
d-wave spin singlet superconducting state on the basis of 
the BdG equations. It has been shown that the AFM or- 
der coexists with the FFLO superconducting state even 
when the AFM order occurs neither in the normal state 
nor in the BCS state. In other words, the AFM phase 
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can be confined in the FFLO phase at high fields, con- 
sistent with the experimental results for CeCoIns. 25,26 
Magnetic instability is enhanced in the inhomogeneous 
Larkin-Ovchinnikov state because the 7r-phase shift of 
the pairing field gives rise to the Andreev bound states 
and produces the large local DOS at zero energy. The 
mixing with the 7r-triplet pairing also enhances the AFM 
order in the FFLO state. 

The structures of the AFM order, namely, the direc- 
tions of the AFM magnetic moment Maf ( i ) and in- 
commensurability 5 ic , are consistent with the recent 
neutron scattering measurement for the HFSC phase of 
CeCoIns. 26 Both Maf(0 and S\c are perpendicular to 
the magnetic field. It has been shown that the ampli- 
tude of the incommensurability |<5ic| is independent of 
the density of spatial nodes of the modulated supercon- 
ducting order parameter. 30 This is also consistent with 
the experimental result. 26 These results indicate that 
the HFSC phase of CeCoIn 5 is the AFM-FFLO state 
in which the incommensurate AFM order coexists with 
the FFLO superconductivity. To obtain unambiguous ev- 
idence for this scenario, we propose another neutron scat- 
tering experiment that can detect the spontaneous trans- 
lation symmetry breaking along the magnetic field. 
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Fig. 7. (Color online) Schematic figures of the H-T-phase dia- 
gram in CeCoIns (a) at ambient pressure (P = Pq), (b) below 
the critical pressure (Po < P < P c ), and (c) above the critical 
pressure (P c < P). 



Observations under pressure may be a further way to 
explore the HFSC phase of CeCoIn 5 . Since the AFM or- 
der is expected to be suppressed by the pressure as in the 
other Ce-based heavy fermions, 49 the AFM-FFLO state 
is gradually suppressed by the pressure, as schematically 
shown in Fig. 7. This is in sharp contrast to the pure 
FFLO state that is enhanced by the pressure as it goes 
away from the quantum critical point. 36 Thus, the AFM 
order should be distinguished by the BCS-FFLO transi- 
tion under pressure below P < P c , as shown in Fig. 7(b). 
When the pressure exceeds the critical value P c , the 
AFM-FFLO phase vanishes, as shown in Fig. 7(c). The 
experimental data for the pressure dependence are con- 
sistent with the enhanced FFLO phase, 19 but the mag- 
netic order in the low-pressure region has not yet been 
studied. 

Finally, we discuss the experimental results for 
CeCoIns. One of the key experiments is that on the 
pressure dependence of the iJ-T-phase diagram 19 men- 
tioned above. It has been shown that the HFSC phase of 
CeCoIns is enhanced by the pressure. This experimental 
result is hardly understood by regarding the HFSC phase 



as an AFM ordered state in the uniform superconducting 
state or in the simple Abrikosov vortex state. Therefore, 
the other quantum condensed state likely emerges in the 
HFSC phase of CeCoIn 5 , such as the FFLO superconduc- 
tivity investigated in this paper or the 7r-triplet pairing 
(pair density wave) proposed by other authors. 31,32,34 
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